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Abstract 

We present a model for next-to-leading order resummed threshold form factors based on a time-like coupling 
recently introduced in the framework of small x physics. Improved expressions for the form factors in A^-space 
are obtained which are not plagued by Landau-pole singularities, as the included absorptive effects - usually 
neglected — act as regulators. The physical reason is that, because of faster decay of gluon jets, there is not 
enough resolution time to observe the Landau pole. Our form factors reduce to the standard ones when the 
absorptive parts related to the coupling are neglected. The inverse transform from iV-space to x-space can be 
done directly without any prescription and we obtain analytical expressions for the form factors, which are well 
defined in all x-space. 
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1 Introduction 



Resummation of large infrared logarithms in form factors and shape variables is essential in order to predict 
accurate cross sections in many phenomenologically relevant processes El El El El ■ In this paper we present 
a model for next-to-leading order (NLO) resummed form factors based on the time-like coupling recently intro- 
duced by B. Ermolacv, M. Greco and S. Troyan in 6 in the framework of small x physics (see also jTj). The 
usual expression for resummed threshold form factors in A^-space is EH^E]: 



/at (as) = exp / 

"'0 



dz- 



Bi asiQ^il -z)) + ■■■ + Di asiQ^l - z)^) + •■•}. 



(1) 



Q is the hard scale of the process. Ai, A2, ■ ■ ■ Bi, ■ ■ ■ Di, ■ ■ ■ are the first coefficients of the functions A (as), 
B (as) and D (as): 



A (as) = Aias + A2a 



B (as) = -Bias + B2C 



D (as) = Dias + 



(2) 



A (as) describes the emission of partons which are both soft and coUinear, B (as) describes hard and coUinear 
partons while D (as) describe partons which are soft and at large angles. The knowledge of the quantities Ai, 
A2, Bi and Di is needed for resummation at next-to-leading order (see later for definition). For instance, in 
the case of the thrust distribution we have |5]: 



Ai - 



2Cf 



A. = 



2Cf 



^ ,67 

' 36 ~ 12 



(3) 



where Ca = -^c, Cp = {Nq — l)/(2iVc), = 1/2, Nc = 3 is the number of colors and np is the number of 
active quark flavors'^. For heavy flavor decays we have instead [5lllUj: 



^1 = ^, 

TT 



A2^^ 



^ ,67 

' 36 ~ 12 



—npTn 
9 



TT 



(4) 



Generally, one assumes as ^ 1 and uses truncated as expansions for the functions A(as), B{as) and D{as) 
in eq. Since the running coupling as (fc^) is integrated over all gluon transverse momenta kt from the 
hard scale Q down to zero, the Landau pole is hit. Therefore as diverges inside the integration region — it is 
certainly not small there — indicating that a truncated expansion for the A, B and D functions is not correct. 
We have conceptually a breakdown of the scheme: resummed perturbation theory assumes as small at the 
beginning, while it ends up with a large as- A prescription has to be assigned to give a meaning to the formal 
expression (Q. 



The standard solution is to expand the exponent of eq. J^l in a function series of the form |3]: 



n=0 



A = /3o as(Q') L, 



exp [L gi (A) -f- 32 (A) -t- as ffa (A) 



L = logA^ 



/at (as) = exp L gi (A) + ^ ag g„+2 (A) 

where 

and /3o = (11/3 Nc ~ 2/3 np)/{A:-K). The functions gi (A) have a power-series expansion: 

00 

9^ (A) = 

n=l 

The effects of the Landau pole in this framework are the following: 

^The value of A2 is given in the MS scheme for the coupling constant. 



(5) 



(6) 



(7) 
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1. the series in eq. (O is divergent as the higher order functions have factorially growing coefficients |1UII11| : 

2. the functions gi{X) have branch cuts starting at A = 1/2 and going up to infinity. The form factors are 
then formally well defined up to a critical value N^rit ~ Q/-^ (A is the QCD scale), above which they 
acquire a (completely unphysical) imaginary part. 

The first problem is solved by means of a truncation of the series to its first few terms — typically two or three — 
which is effectively a prescription of the Landau pole. This is the so-called fixed logarithmic accuracy, which we 
will describe later. The resummation to leading order requires the knowledge of the function gi, the resummation 
to next-to-leading order also requires the knowledge of g2, and so on. As far as the second problem is concerned, 
one simply restricts himself to a fiducial region in A^-space below Ncrit- The inverse Mellin transform from A^- 
space to x-space, at next-to-leading log(l — x) accuracy, is well defined up to Xcrit ~ 1 — A/Q, above which 
singularities occur. However, a form factor formally well defined in the whole a;-space and containing all the 
requested log(l — x) terms can be obtained by means of an additional prescription for the contour integration 
in A^-space, the so-called minimal-prescription jl2) . 

Another common solution to this problem involves the renormalon calculus The latter uses the in- 

consistency discussed above to get information about non-perturbative effects, in the form of power-suppressed 
corrections to the cross sections: ^ 

6a^c(^) . (8) 



The exponent b of the power correction can be computed, but not its coefficient c. The resulting information 
is therefore of a qualitative kind and, ultimately, perturbation theory is deprived of predictive power. In fact, 
one can substantially modify the spectra coming out of a renormalon calculation by changing the prescriptions 
in the Borel plane, which are completely arbitrary. 

Our perspective is different: we want to use resummed perturbation theory in a predictive way, by curing 
the " disease" of usual resummation formula. Therefore we begin by re-analyzing the derivation of the standard 
formula. as{k^) occurs in eq. (Q) because we are computing a so-called inclusive-gluon-decay quantity, in 
which one does not observe the development of the jets originating from the gluons emitted by the hard 
partons. One then sums over all possible final states, i.e. formally over all cuts of dressed gluon propagators, 
reconstructing their discontinuity. The radiated gluons have a positive virtuality fc^ > 0, as the development of 
jets is intrinsically a time-like process. Absorptive effects, i.e. the —iir terms in gluon polarization functions, 
are usually neglected in literature: by taking them into account, one enforces that all gluons participating to 
the cascade are unstable particles [HI- They are analogous to the quasi-particles of statistical physics, possessing 
in lowest order — non interacting quasi-particles — a non-zero width . The outcome of the inclusion of this 
physical effect is that as(fc^) does not occur any longer in the resummation formula, and it is replaced by a 
different function of fc^, 

asikl) ^ 55 (fc?), (9) 

which does not possess the Landau pole. We call this new function "effective coupling", as it specifies the 
effective strength of the interaction in the gluon cascade — not in a general QCD process. The main point is 
that the effective coupling never becomes large, as the typical expansion parameter is 

< < 1 for any fc? (10) 

TT /ioTT 

and one never leaves the perturbative domain. By including absorptive effects related to the coupling constant 
— usually neglected — we derive an improved expression for the resummation formula of the form factors. 
The resummation formula is free of Landau pole pathologies, it does not involve any new free parameter and is 
strictly predictive. 
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Figure 1: Coupling constant as a function of gluon transverse momentum squared k. 
coupling as{kf), eq. ((T5|) : dashed line: standard coupling as{kf), eq. ifTHl . 



Solid line: effective 



2 The improved resummation formula 

By taking into account renormalization effects in the form factors, the tree-level, momentum-independent cou- 
pling is replaced by the integral on the discontinuity of the gluon propagator ^1 : 



as 



asik 



1 
2Td 



ds Disc, 



1 



s Po log 



(11) 



The discontinuity is defined as usual: DisCsF(s) = F(s -I- ie) — F(s — ie), with e being a positive infinitesimal 
number. We now close the integration contour with a circle of radius k^ and a circle of infinitesimal radius. 
By using Cauchy's theorem and neglecting the residue of the pole in s = — in order the preserve asymptotic 
freedom ^, we obtain: 

dip 1 



as{kt 



2tt 



Po log j^+if 



(12) 



The standard resummation formula ^ is obtained by neglecting the imaginary term itp in the denominator. In 
hard processes Q 3> A and one expects on physical grounds that kt ^ Q, implying: 



log^ » ^- 

As a consequence, one obtains the coupling evaluated at the gluon transverse momentum squared: 

1 



asikt) 



(3o log ■ 



asik^). 



(13) 



(14) 



However, the assumption H13() is not correct, because the transverse momentum kt is integrated from Q down to 
very small scales. We propose to modify the resummation formula according to this criticism: approximation 
H13|l disregards absorptive effects related to the gluon jet decay, which are important in the infrared region, and 
therefore has to be avoided. By performing the integration exactly, we obtain: 



as{kt) = — 
Po 



11 log ■ 
arctan ■ 

2 TT TT 



(15) 



*The integral on the discontinuity of the gluon propagator is transformed into the integral over a circle of radius plus the 
residue of the pole in s = — A^. In our pragmatic approach, we have just omitted the latter. 
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Figure 2: Form factor in A^-space /at (as) for the beauty case. Solid line: our model; dashed line: standard form 
factor, i.e. r = in our model. 



The effective coupling H15|l approaches the standard one in the asymptotic region fc| A^, but it does not 
contain the infrared pole in fc^ = ^. Furthermore, as{kf) is positive definite, monotonically decreasing in all 
the kt range and has a finite limit at zero momentum (see fig. 1): 



.1™ as(fct) 



1 



(16) 



According to eq. H15|) . the effective coupling deviates from the standard one and saturates at a scale of order 

kt - Ae'^/^ _ ;LGeV, (17) 
for A 200 McV. Our improved expression for the form factor at one-loop approximation reads: 
' _ 1 



/7v(as) = exp / dz — 
Jo 



^3 (1^-) dkl 

Q2(l-z)2 fct 



Ai asik^t) + Bi a5(g'(l - z)) + Di AsiQ^^ " zf)}, (18) 



where the expression for the effective coupling as is given by eq. (|15|l . By comparing with the standard 
resummation formula we see that the only difference is the appearance of as in place of as. Let us now 
consider two-loop effects. Note that eq. (Illll can be written as: 



2ni 



ds DisCj. — — — . 



where 



as(Ai ) 



2\1L 



1 



/3olog^ 



is the one-loop coupling. We now assume that eq. H19|) generalizes to the two-loop coupling 

1 /3i loglog^VA^ 



as (/i 



/3ologMVA2 f3l logVVA^ 



(19) 
(20) 

(21) 



^We define the arctan function as the one being continuous in zero and discontinuous at infinity, with image in the interval 
(-7r/2,+7r/2). 
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Figure 3: /at (as) for the top case. Solid line: our model; dashed line: standard form factor, i.e. r = in our 
model. 

where Pi = (51 — 19/3 ni?)/(87r^). Therefore, the subleading contribution to the effective coupling as{k1) is: 



/3i 



^. 1 -/3iloglog(-s/A2) 
as DisCc 



s PI log2(-s/A2) 



1 1 logfc2/A2"l 

h - arctan 

2 TT TT 



k"^ 1 
log ^ + 2 



TT^ + log- 



A2 



(22) 



It is not easy to derive the terms proportional to A2. Let us make the following guess: 



4(fc?) 



/ dsDisCs ^S-!^ 

1m Jo s 



IL 



1 

2^ 



ds Disc, 



s/32log"(-s/A2) 



/32(^2+log2 ^2/^2)- 



(23) 



An explicit verification of the validity of eq. H23|l would require a three loops computation, which is beyond 
the purpose of the present work. That is because A2 first occurs at order a| and running coupling effects 
involve at least one additional power of as, i.e. a|. in total. Let us stress that also the standard replacement 
a| — > as(fct)^ has never been explicitly checked, as its verification requires a three loop computation as well. In 
eq. (j23|l we just assumed that the s-discontinuity had to be taken after squaring the gluon polarization function; 
it is remarkable that our choice gives a simple analytic result. Other choices are possible at this level, such as 
squaring the effective coupling: a|(fc|) asikf)"^ ^ . Anyway, since two-loop terms are rather small, our choice 
seems acceptable in a model of soft gluon dynamics. 

In conclusion, our expression for the rcsummcd form factor at next-to-leading order reads: 
/Ar(as) = exp / dz— / -^[A^asikf) + A2al{kf) + ■ ■ ■] + 

Jo '-Jq2(1_^)2 fcj 

+ Bi5s(Q'(l -z)) + -- - + D^asiQ^l - zf) + •••}, (24) 

®With this choice there is only one effective coupling, while with our choice there are two effective couplings. 
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where: 



1 

Pi 



1 1 ^ log^ 
arctan — 

2 TT TT 



1 



PS + log^ fc2/A2 



1 logfc2/A2 
— arctan 

TT TT 



log 



A2 



■log 



log 



A2 



+ 1 (25) 



and is given by eq. (|^ . 



Eq. (|24|l is the main result of this work and replaces eq. The main problem of eq. discussed before, 
does not occur anymore in our improved expression: the effective couplings as and a| do not diverge for small 
values of the argument and actually remain acceptably small. 



Our treatment of the running coupling constant resembles the exponentiation of the tt^ terms in the K- 
factor of Drell-Yan coefficient function It is well known that, going from the space-like kinematics of 

Deep- Inelastic-Scattering {q^ < 0) to the time- like kinematics of Drell-Yan (g^ > Q), double logarithms in the 
vertex correction generate a large constant term proportional to tt^ , because: 



Since double logs exponentiate, it has been suggested that also the tt^ terms exponentiate: 

g-cas log^(-(J^-iO) ^ gC as TT^ g-c as log^(i}^) 



(26) 



(27) 



where c is a positive constant ^. It has been however established that the exponentiation of tt^ is not exact, 
being violated at two loops by terms of the form C(2)log^(— — jQ), where C. is the Rieman zeta function 
(C(2) = TT^/G) PZj. Analogously, in our case, we cannot expect to have a complete resummation; we resum, 
however, a set of constants which are certainly there, and this seems to be a reasonable approximation. 



3 The functions gi and g2 

Let us begin this section by defining the logarithmic accuracy in our model. We have powers of the coupling at 
the hard scale, 

a = asiQ^), (28) 



multiplied by powers of the infrared logarithm, L = log A^, i.e. monomials of the form 



a L 



(29) 



Within the usual resummation scheme. 



• by leading logarithms {LL) we mean all the terms in the exponent of the form factor of the form 

a" L"+i with n = 0,l,2,---oo; (30) 

• by next-to-leading logarithms (NL) we mean all the terms of the form 

a" with n = 0,l,2,---cx). (31) 



our work we are dealing with a geometric series (cfr. eq. 1121 *). while in the /C-factor an exponential series is resummed. 
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The coupling as{kf) at one-loop, when expanded in powers of a, produces a series of leading logarithms: 

asikf) = a - /3o«2 log ^ + p',a' log^ ^ + ' ' ' ^^^^ 

The powers of log /Q'^ give indeed powers of L = logiV after the integration over kt and z in eq. 
According to eq. 115|) . the relation between the standard and the effective coupling reads ^: 

as(kt) ^ ^T- arctan\TT/3Qas(kf)] for as(k?) > 0. (33) 
7r/3o 

Therefore, the expansion of our effective coupling in the parameter a is: 

as{k^) = a- Poa' log ^ + P^a' log' ^ - ^-^a' + ■■■ (34) 

Comparing eq. with eq. we see that the expansion of the effective coupling involves subleading 

logarithms compensated by powers of /Jqtt, coming from the absorptive effects discussed before. The main point 
is that the Pqtt terms have a fundamental regulating effect over the logarithmic ones; therefore one has to 
consider them on equal footing in the counting. The hierarchy in our model can therefore by defined as follows: 



• by LL we mean all the terms of the form 

a" iPoL)"+^~-'' (Pott)'" with fc = 0, 1, ■ • • n + 1; (35) 

• by NL we mean all the terms of the form 

a'' iPoLY'-'' ipon)'' with k = 0,l,---n. (36) 



The usual scheme is "minimal" in the sense that it deals only with powers of the logarithms and of the coupling; 
the related problem is that of the Landau pole discussed above. Instead, our scheme, in order to solve the 
Landau problem, has to be non-minimal. 

The next step is to perform the integrations within the above defined next-to-leading accuracy. In general, we 
believe that fixed logarithmic accuracy is a consistent approximation in our scheme because series of logarithms 
of higher order are suppressed by powers of as, which is always small. The integrations over transverse and 
longitudinal gluon momenta are easily done within the standard approximation 4 : z^~^ — l ~ —9 {1 — z — 1/n) , 
where 9 is the step function and n = Ne^'^ , with being the Euler constant. This approximation, as shown 
in our previous paper 10 , misses for instance terms of the form: 7r^a"i"^^, which are leading according to 
our power counting. However, in our model, we are interested in picking up only those terms proportional to 
/3o7r, which come from the analytic continuation of the coupling from space-like to time-like region. These are 
the terms coming from the effective coupling, having the regulating effect on the Landau pole. 

In order to obtain the functions gi and g2, we express the logarithm of the hard scale logQ^/A^, coming 
from the integrations, in terms of the coupling evaluated at a general renormalization scale /i by means of the 
formulas: 

and 

®One has to use the relation arctan(l/x) = 7r/2 — arctan(a;), which is valid for a; > with our definition of the arctan function. 
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We finally Taylor expand the exponent on the r.h.s of eq. H18|l in powers of 7^;, (3i and log/i^/Q^ up to first 
order included. Our result for the leading function gi reads: 



Ai 



- (1 - 2 A) log[(l -2Xf+ r^] + 2 (1 - A) log[(l - A)' + r^] - log[l + r^] + (39) 



-TT A^ + (1 - r^) arctan - + (1 - 2A) - 
r I 



1 -2A „ 
arctan 2 



(1 - A)^ - 



arctan ■ 



1 - A 



The subleading function 52 reads; 



32(A;r) 



^(log[(l-2A)2+r2]-log[l + r2] + - 
4/3o I r 



Bi 
2/3o 

A2_ 

2 



log[(l-A)2+r2]-log[l 



1 



1 1 — 2A 

-TT A + arctan (1 — 2A) arctan 

r r 

1 , , 1 - A 

— TT A + 2 arctan 2 (1 — A) arctan 

r r 



2 ^ log[(l - 2A)2 + - 21og[(l - \f + + log[l + r^] + 



r 

16 /33 



1-2A , 1-A 1 

— (1 — 2A) arctan V 2(1 — A) arctan arctan - 

r r r 



i|- log2[(l - 2A)2 + r^] + 21og2[(l - Xf + r^] - log2[l 

3 I 



+ ^(1-2A) 
r 



2 1 - 2A 
- 1 H — arctan 



2 1-A 

1 arctan 

TT r 



log[(l-2A)2+r2] + 



27r 



log[(l - A)^ + r"] + — 
r 



2 1' 
-1 H — arctan — 



log[l + r2]^ 



4 arctan - 



1 - 2A 



-TT + arctan - 



1 - 2A 



8 arctan ■ 



1-A 



TT — arctan ■ 



1 - A 



4 arctan - 



-TT + arctan ■ 



^ log ^ log[(l - 2A)2 +r']-2 log[(l - A)^ + r^] + log[l + r^] + 



2 

r 

Ai7e 
2/3o 



1 - 2A , , 1-A 1 
— (1 — 2 A) arctan h 2(1 — A) arctan arctan - 



r r 
^2 , 21 i„„r/i \\2 , 21 



log[(l - 2A)^ + H] - log[(l - A)^ + H] + 



1-2A , , 1-A ttA 

— (1 — 2 A) arctan h (1 — A) arctan — 

r r 2 



(40) 



We have defined — note the general renormalization scale fi in the coupling: 

A = asifJ-'^) L and r = tt /?o Q;s(/i^). 



(41) 



Note that A and r are obtained from each other by interchanging L with tt. All this is in line with the logarithmic 
accuracy defined above. 

A few remarks are in order. The first one is that, unlike the standard case, gi and 32 do not depend only on 
A but also on r, i.e. on as alone. Our gi and 32 are then "non-minimal" , as a consequence of the inclusion of 
absorptive constants in the coupling. The parameter r 7^ acts as a regulator of Landau-pole singularities: gi 
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Figure 4: Renormalization scale dependence of form factor in A^-space /at (as) for the beauty case. Solid line: 



Q^; dashed line: /i^ = dotted hne: fj.'^ = iQ^. 



and g2 are regular for any value of A. As expected, in the limit r — > we recover the standard gi and g2 0]: 

Ai 



and 



ffi(A;r = 0) 



2/3oA 



{-(1-2 A) log[l-2A] + 2 (1-A) log[l-A]} 



(42) 



52(A;r = 0) = 



Di log(l - 2 A) Bi log(l - A) A2 [log(l - 2 A) - 2 log(l - A)] 

2/3o Po 2/30^ 

Ai Pi [2 log(l - 2 A) + log^(l - 2 A) - 4 log(l - A) - 2 log^(l - A)] 

Ai log(/iVQ') [log(l - 2 A) - 2 log(l - A)] , Ai je [log(l - 2 A) - log(l - A)] 



(43) 



2/3o 



If we expand our 171 and 172 for small r, we find that the leading corrections are of order r^A": this means cor- 
rections of next-to-next-to-leading order (NNLO) a^'^^L'^ in the standard counting, which had been overlooked 
in our previous evaluation of the 173 Our gi and 172 contain some terms from the standard 53, 54, 35, ... or 
from the coefficient function C{as) multiplying the form factor. It is remarkable that eqs. (|39|l and (|4()|l involve 
only powers of the log and arctan functions. 

In fig. 2 we have plotted the form factor /at (as) as a function of N for the decay of a beauty quark, i.e. 
for Q = rrib^, for which as = 0.21 and np — 3. The dashed line is the plot of the standard form factor, i.e. of 
the same form factor in the limit r ^ 0. We have plotted moments up to iV = 20, because above this point 
the standard form factor, as already discussed, becomes singular. Our form factor is reasonably close to the 
standard one up to A*" ~ 10, the latter being steeper for larger A^. For iV = 20 there is a difference of a factor 2 
circa. Fig. 3 shows similar plots for the top case, i.e. for Q — rrit for which as ~ 0.11 and np — 5. Due to the 
increased hard scale by more than one order of magnitude: 



• differences between the two models are much smaller; 

• the peak above one, related to the subleading, single-logarithmic effects, is barely visible, while it is rather 
pronounced in the beauty case. 

^In general, for heavy flavor decays Q = 2Ex, where Ex is the energy of the hadronic final state. One can set 2Ex = "ib in 
b — > s-y, while this is not true in semi-leptonic b ^ u decays 1181 . 
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f (x;as) 

5- 
0- 
5- 



0.8 0.85 0.9 0.95 

Figure 5: Form factor in a;-space /(x; as) for the beauty case. Solid line: our model; dashed line: standard form 
factor. 

f (x;as) 




X 



0.98 0.985 0.99 0.995 



Figure 6: Form factor in a;-space f{x;as) for the top case. Solid line: our model; dashed line: standard form 
factor. 

The scale dependence is shown in fig. 4, where we have plotted /Ar(as) for the beauty case for fi^ = Q^, 
fj? — and /i^ = AQ"^. The scale ambiguity is smaller than the difference between the models shown in 

fig. 2. 

Let us now discuss the extension of our model to NNLO, i.e. the computation of the function 173. It is 
clear that one has to make similar guesses to the ones needed at NLO to evaluate the terms proportional to 
A2. Apart from that, the computation, although technically rather cumbersome, does not seem to present any 
specific difficulty. In general, our fixed logarithmic accuracy for the exponent in the form factor, 

J: ^ L gi{X;r) + g2{X;r) + a g3{X;r) + a^g4{X-r) + ■ ■ ■ , (44) 

is an expansion with better convergence properties than the usual one, because the effective coupling becomes 
at the most of order one. In other words, our expansion is certainly better than an expansion in l/(/?o7r)", 
which is already convergent as l/(/?o7r) < 1. We expect the extension of our model to NNLO to be relevant 
also for beauty physics. In a previous work ^U] we have found instead that usual NNLO effects could not be 
included in the case of beauty decay. In fact, the hard scale Q — mi, was not large enough to avoid effects of 
the divergence of the perturbative series even at relatively low values of x. Such divergence was related to the 
integration over the Landau pole and is therefore absent in the present scheme. 
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Figure 7: Renormalization scale dependence of form factor in x-space f(x; as) for the beauty case. Solid line: 
/x^ = Q^; dashed line: = (3^/4 ; dotted line: /i^ = 4(3^. 

4 Form Factors in x space 

Up to now we have considered form factors in iV-space. The A^- moments are indeed physical quantities, but in 
practice a measure of the moments for large N is difficult. Therefore, let us transform back to x-space. The 
inverse transform of /at (as) is defined as: 

f(x-as)^— dNx~^fN{as), (45) 

where C is a constant such that the integration contour lies to the right of all the singularities of /jy. A standard 
computation gives to next-to- leading log(l — x) accuracy jj]: 

-S) ^-^Ul-x- o) exPKg:(r;0+,.(r.)] | 
\ r[l-^(rgi(r;r))] J 

where we have defined: 

T = l3o asin'^) I and I = -log(l-j:). (47) 

In eq. H46|l we have neglected terms 0(1 — a;), which are small in the threshold region. The form factor f{x; as) 
in eq. is plotted in fig. 5 for the beauty case. Also shown is the standard form factor, i.e. the same formula 
for r = 0. Our form factor is formally well defined up to x = 1, while the standard one presents a singularity at 
Xcrit ~ 1 — A/Q, SO it is plotted only below this point. There are substantial differences between the models: 
our distribution is broader, the peak is smaller and occurs at much larger x than in the standard case. The 
reason is the following. Our model contains an effective coupling which approaches a constant value for iV — > oo 
or, equivalently, x — > 1 (see eq. (I16|| V Therefore, for x ~ 1 our model resembles a frozen coupling scheme, while 
the standard form factor contain a coupling divergent at the non-zero momentum A — the well-known Landau 
pole. In fig. 6 we plot /(x; as) for the top case in both models; as expected, differences are much smaller with 
respect to the beauty case. The dependence of our /(x; as) from the renormalization scale /x is shown in fig. 7 
for the beauty case and turns out to be moderate. 



5 Phenomenology 

Let us now discuss the phenomenological relevance of our model. A natural application concerns beauty physics, 
i.e. a case with a moderately large hard scale — when the hard scale is very large, as in top decays, there are 
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basically no differences between our model and the standard one. In the beauty case, large non-perturbative 
corrections are expected. The main problem, in our opinion, is that of combining perturbative and non- 
perturbative effects in the most efficient way. In many phenomenological studies, one assumes a shape for the 
non-perturbative effects and convolutes it with a frozen-coupling distribution. Since asijub) ^ 0.21 (that is, 
rather small) and the coupling docs not increase at small kt, the perturbative effects turn out to be extremely 
small. Since the theoretical prediction is a convolution of perturbative and non-perturbative effects, it is 
completely dominated by the assumed form of the latter. The standard rcsummation scheme — running 
coupling — includes much more perturbative dynamics, but it is technically more difficult to implement in a 
phenomenological analysis and it has been barely used. That is because the perturbative form factors in x-space 
are generally computed with a numerical contour integration, by using the so-called minimal prescription jl2| . 
This is a consistent prescription based again on " minimality" : the correct logarithms of 1 — a: are included and 
the distribution is formally well defined up to a; = 1, but no physically motivated effects are included, only 
logarithms. Our approach includes instead a (perturbative) physical effect, namely parton decay, giving rise 
to analytic distributions well-defined in the whole x-space. One can then easily convolve our spectra with an 
assumed form of non-perturbative effects: 



where fnp is a non-perturbative function and f^x is the spetrum to be compared with the data. Cfe denotes 
schematically non-perturbative parameters entering fnp. 

It is well known that non-perturbative effects broaden the peaks of perturbative spectra and shift them 
to smaller x's. They have, as intuitively expected, a smearing effect, representing the final stage of parton 
evolution, not described by perturbation theory. We have not done yet any detailed phenomenological analysis, 
but the following points can already be made. Our model gives a broader spectrum compared to the usual one, 
while it gives a peak at larger values of x. The peak we find corresponds indeed to a jet invariant mass below 
1 GeV. We then expect that our model has to be supplemented with non-perturbative effects shifting the peak 
towards smaller x. 



6 Conclusions 

Our main result is the improved threshold resummation formula given by eq. (|24|l . The latter involves — in 
place of the usual running coupling possessing the Landau pole — an effective coupling which is not singular in 
the infrared region and remains rather small in all the integration domain. It approaches the standard coupling 
for large transverse momenta, while it deviates from the latter at a scale of order 



The effective coupling, unlike the usual one, includes absorptive effects related to the decay of gluon jets. The 
physical explanation of the absence of the Landau pole in our resummation formula is the following. The 
inclusive gluon branching, inducing the running coupling in resummation formula, is described by taking into 
account absorptive effects in addition to the usual dispersive ones. The gluon cascade is then treated, more 
realistically, as an intrinsically unstable process and there is not enough resolution time to see the effects of the 
Landau pole. 

By performing the integrations over longitudinal and transverse momenta with next-to-leading accuracy, we 
have obtained the functions gi and g2 in the exponent of the form factor given in eqs. H39|) and 14U|) respectively. 
These functions define next-to-leading order resummed threshold distributions which are well defined in the 




(48) 



Ae""/^ - 1 GeV. 



(49) 
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whole moment or physical space. There are large differences between our next-to-leading order form factor and 
the standard one for beauty decays. We believe that our model can be used to describe the variety of spectra 
in inclusive or rare b decays. Natural extensions of our approach involve resummation of shape variables or the 
inclusion of mass effects in semi-leptonic b ^ c transitions. 
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